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SUMMARY 


A structural performance and resizing (SPAR) finite-element computer program and 
NASA structural analysis (NASTRAN) finite-element computer programs were used in the 
thermal stress analysis of the space shuttle orbiter subjected to reentry aerodynamic 
heating. A SPAR structural model was set up for the entire left wing of the orbiter, 
and NASTRAN structural models were set up for (1) a wing segment located at midspan 
of the orbiter left wing and (2) a fuselage segment located at midfuselage. The 
thermal stress distributions in the orbiter structure were obtained and the critical 
high thermal stress regions were identified. It was found that the thermal stresses 
induced in the orbiter structure during reentry were relatively low. The thermal 
stress predictions from the whole wing model were considered to be more accurate 
than those from the wing segment model because the former accounts for temperature 
and stress effects throughout the entire wing. 


INTRODUCTION 


The space shuttle orbiter is designed to be flown as many as 100 missions with- 
out excess thermal and mechanical loadings. To establish confidence in the orbiter 
thermal and structural integrity, it is essential to fully understand both the 
thermal and structural performance of the orbiter subjected to reentry aerodynamic 
heating and aerodynamic loading. 

Because the number of onboard thermocouples is extremely limited, it is impossi- 
ble to generate accurate temperature distribution within the orbiter structure based 
on thermocouple data for estimation of the temperature in each structural component. 
For this reason analytical thermal analysis (for example, finite-element heat trans- 
fer analysis) of the orbiter is necessary, ohe thermal analysis can give relatively 
accurate temperature distribution in the orbiter structure, making it possible to 
determine the temperature level in each of the orbiter structural components. This 
analysis can show if the design limit temperature of 350°F is exceeded. (Heating 
beyond 350°F will certainly degrade the aluminum structural material.) Addition- 
ally, a thorough knowledge of the structural temperature distribution is necessary 
for an accurate thermal stress analysis. 

The flight load data obtained from onboard strain gage measurements contain both 
the thermal and mechanical stresses. Unfortunately, these two stress components are 
not easily separated experimentally. To obtain the mechanical stresses, the thermal 
stresses must be removed from the strain-gage-measured stresses. This can be done 
analytically by first calculating the thermal stresses and then removing them from 
the strain-gage-measured stresses to give the true mechanical stresses. For the 
thermal stress calculations, the structural temperature distributions obtained from 
the heat transfer analysis may be used as input to a structural model. Extensive 
work on the heat transfer analysis of the orbiter was conducted by Ko, £iinn, and 
Gong (refs. 1 to 7). 

The purpose of this report is to use the finite-element method to calculate 
thermal stresses in the orbiter structure using the structural temperature distribu- 
tions obtained from references 4 to 7 as thermal loadings to the structural models. 
Finite-element structural models were set up for the entire left wing, a wing segment 



located at midspan of the left wing, and a fuselage segment located at midfuselage. 
Thermal stress distributions in the orbiter structure were obtained, and the crit ’ 
cal high-stress regions were identified. These analyses also provide a baseline 
establishing element mesh sizes which will be adequate for thermal stress analysis 
of large aerospace structures - 
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NOMENCLATURE 

quadrilateral membrane and bending element 
two-node tension/compression/torsion element 
bar element 

zero-length element for elastically connected geometrically 
coincident joints 

triangular membrane element 

quadrilateral membrane element 

quadrilateral shear panel element 

flexible felt reusable surface insulation 

fuselage segment structural model at station Xq877 

high-temperature reusable surface insulation 
joint location (or node) 

low-temperature reusable surface insulation 

NASA structural analysis 

room temperature vulcanized 

strain isolation pad 

structural performance and resizing 

space transportation system, flight 5 

thermal protection system 

whole wing structural or thermal model 

wing segment structural or thermal model at wing station | Y 0 j 240 
station in x direction 
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z 0 

x, y, z 
Aa 



station in y direction 

station in z direction 
rectangular Cartesian coordinates 
change in angle of attack, deg 

normal stress in x direction, lb/in z 
normal stress in y direction, lb/in^ 
shear stress, lb/in^ 


DESCRIPTION OF PROBLEM 


Figure 1 shows a composite top and bottom view of the space shuttle orbiter. The 
orbiter parts selected for the present study are the major load-carrying regions of 
the left wing, and the midfuselage cross section (fig* 1)* The problem i 3 to calcu- 
late thermal stresses induced in the space shuttle orbiter structure during reentry 
aerodynamic heating. The thermal loadings to the orbiter structure were based on the 
space transportation system flight 5 (STS-5) reentry trajectory shown in figure 2 
(ref. 4). Three regions of the orbiter were selected for the thermal stress analy- 
sis. The first region was the whole left wing bounded by chordwise stations Xo1040 
and Xol 365. The second region was a segment of the left wing bounded by spanwise 
stations Yo~240 and Yq- 254. The third region was a fuselage segment bounded by fuse- 
lage stations Xo877 and X()880. The elevon, leading edge regions, wheel well door, 

and bay door were not included in the analysis because the present study concerns 
only the major load-carrying structures. Also the thermal protection system (TPS) 
was excluded in the thermal stress analysis because it is not a load-carrying struc- 
tural component. However, for the heat transfer analysis, the vheel well door, the 
landing gear, the bay door, and the TPS had to be included. 


DESCRIPTION OF STRUCTURES 


Wing 

As shown in figure 3, the ribs of the wing are aluminum truss systems made of 
vertical, horizontal, and diagonal members pin- joined together. All of the spar 
webs, the wheel well vertical walls, and the wing root vertical wall (fuselage wall) 
are ma4e of corrugated aluminum panels. Both the lower and upper wing skins lying 
between stations Xol 191 and Xol 365 are made of hat-stringer-reinforced aluminum 
panels. The lower and upper wing skins lying between stations Xol 040 and Xol 191, 

and the leading edge beam assembly are made of aluminum honeycomb-core sandwich 
panels. The landing gear door is made of double-walled, hat-stringer-reinforced 
aluminum panels separated by aluminum stringers. The entire lower wing surface 
is covered with high-temperature reusable surface insulation (HRSI) tiles, with a 
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layer of strain isolation pad (SIP) lying between the wing skin and the HRSI for the 
absorption of strain incompatibility between the aluminum wing skin and the HRSI. 
Most of the upper skin near the leading edge region is covered with low- temperature 
reusable surface insulation (LRSI) tiles. A SIP layer lies under the LRSI to func- 
tion the same as the SIP does for the HRSI. The rest of the upper wing skin, which 
is subjected to low heating, is covered with highly flexible felt reusable surface 
insulation (FRSI), under which there is no SIP layer. Room temperature vulcan- 
ized (RTV) rubber bonding agent was used in bonding the thermal protection system 
(TPS) to the wing surfaces. Some of the gaps between the TPS tiles in the high- 
temperature regions are filled with ceramic-coated aluminum mat (gap fillers) to 
prevent hot gases from coming in contact with the substructure at the bottom of 
each gap. Figure 4 shows the structural details of the wing segment of the 
orbiter right wing. 


Fuselage 

Figure 5 shows the fuselage cross section located at station Xq 877. Both the 

fuselage bottom- (or belly) and the sidewall are made of T-stif fener-reinforced 
aluminum skins. The lower and the upper glove skins (except for the leading edge 
region) are made of hat-stringer-reinforced aluminum skins. The leading edge region 
of the glove skin is an aluminum honeycomb-core sandwich structure. The bay door is 
a sandwich structure made of honeycomb core and graphite-epoxy composite skins. A 
small portion of the bay door inner surface is covered with a layer of RTV rubber to 
serve as a heat sink. The fuselage bottom, lower glove, glove leading edge region, 
and part of the glove upper surface (near the leading edge region) are covered with 
HRSI. Most of the upper glove outer surface is covered with LRSI. The lower por- 
tion of the sidewall outer surface is covered with FRSI, and the upper portion with 
LRSI. The upper outer surface of the payload bay door is covered with a layer of 
FRSI, and lower outer surface with LRSI. 


FINITE-ELEMENT MODELING 


Because of the complex nature of the space shuttle orbiter structure, some 
structural simplifications were necessary before setting up the finite-element 
structural models so that the computations would be manageable. Excessively 
detailed models could lead to excess computation time, for which the gain in solu- 
tion accuracies might not be high enough compared with the solutions obtained from 
simpler yet reasonably detailed models. The previous heat transfer analysis of the 
orbiter (refs. 1 to 7 ) showed that representing the hat-stringer-reinforced skins 
(wing and glove skins), T-stif fener-reinforced fuselage skin, and honeycomb-core 
sandwich skins with smooth panels of effective thicknesses could give sufficiently 
accurate temperature solutions. Therefore, in setting up the structural models for 
the orbiter, a similar approach was adopted. 

For the thermal stress analysis of the whole orbiter left wing, a structural 
performance and resizing (SPAR) finite-element computer program (ref. 8) was used; 
for the thermal stress analysis of the wing and fuselage segments, a NASA structural 
analysis (NASTRAN) computer program (ref. 9) was used. 
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Whole Wing 


The SPAR finite-element structural model (WING) setup for the entire orbiter 
left wing is shown in figure 6. This wing structural model was obtained from direct 
modification of the whole wing thermal model used earlier for the heat transfer anal- 
ysis by Ko, Gong, and Quinn (ref. 4). The wing upper and lower skins were modeled 
with quadrilateral membrane (E41 ) and triangular membrane (E31 ) elements. Only two 
elements were modeled between spar caps, and only one element between the rib caps.- 
The anisotropic material properties were used to account for the effect of the hat- 
stringers. The fuselage wall (or wing root wall) was modeled with E41 elements. 

The spar and rib webs and the loading edge beam assembly were modeled with quadri- 
lateral shear panel (E44) elements. The spar and rib caps and the rib truss members 
were modeled with two-mode bar elements (E23) which have only axial stiffness. The 
wing root nodes are elastically connected to geometrically coincident points through 
zero- length E2S elements to simulate the connection of the wing root to the fuselage 
structure. The size of the entire wing structural model is as follows: 


232 
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elements 

67 

E44 
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Wing Segment 

. The NASTRAN finite-element structural model (WS240) for the orbiter wing segment 
bounded by wing stations Yq- 240 and Yg-254, is shown in figure 7. The upper and 

lower wing skins, spar webs, and rib cap webs werfe modeled with quadrilateral mem- 
brane and bending elements (CQUAD2) • The spar caps and the rib truss members were 
modeled with two-node tension-compression- torsion elements (CROD). 

In order to approximate the actual deformation field of the whole wing, two 
boundary conditions were used for the WS240 structural model. The Yo~254 plane 
was fixed (no displacement in the y direction), but the displacements in the x 
and z directions were permitted. For the Yq- 240 plane, two boundary conditions 
were used: 

1. Plane stress — No rotations with respect to the x, y, and z axes, 
but free to move in the x, y, and z directions. 

2. Plane strain deformation — The y displacements for all the nodes 
lying in the Yo-240 plane were set identically, and rotations with respect 

to the x, y, and z axes were constrained. 

The size of the entire WS240 structural model is as follows: 

204 grid points 

121 C QUAD 2 elements 

1 39 CROD elements 
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Fuselage Segment 


Figure 8 shows the NASTRAN finite-element structural model (FS877) setup for the 
fuselage segment, bounded by the two planes at fuselage stations Xq 877 and Xq880. 

The bay door was omitted because it is not a major load-carrying structure. The 
fuselage segment was modeled with CQUAD2 and CROD elements. The Xq 877 plane was 

fixed (no displacement in the x direction), but the displacements in the y and z 
directions were permitted. The deformation of the Xq 880 plane was constrained to 
be plane strain deformation. The entire FS877 structural model has 

62 grid points 

89 CQUAD2 elements 

9 CROD elements 


THERMAL LOADINGS 


Figure 9 (taken from ref. 4) and figure 10 (taken from ref. 5) respectively show 
the structural temperature time histories for the wing skins and the fuselage skin. 
Notice that the structural temperatures for most of the wing skin and fuselage skin 
stations reached their respective peak values at t = 1700 sec from start of reentry. 
In this analysis, the structural temperature distributions at t = 1700 sec were used 
as thermal load input to the structural models for the thermal stress calculations. 

For both WS240 and WING structural models, the structural temperature distribu- 
tions were obtained from the computer outputs of the earlier heat transfer analysis 
conducted by Ko, Quinn, and Gong (ref. 4). For the FS877 model, the structural tem- 
perature distribution was obtained from the computer outputs of the heat transfer 
analysis done by Ko, Quinn, and Gong (ref. 5). Figure 11 shows the thermal loadings 
at wing station | Yq | 240 for the WS240 and WING structural models. The difference in 

the shape of the calculated structural temperature distributions in figure 11 is 
caused by the number of elements in the two models. Similar distribution can be 
expected in the thermal stress calculations. Figure 12 shows the thermal loading 
used for the FS877 structural model. 


RESULTS 


Figures 13 to 15 respectively show the distributions of the chordwise stress 
o x , spanwise stress Oy, and the shear stress T X y in the orbiter wing lower skin 
calculated from the WING structural model. The peak compression of both a x 

(= -9097 lb/in 2 ) and a y (** -2897 lb/ in 2 ) occurred near the wing root. The peak 

tension of o x = 2836 lb/in 2 and o y * 2405 lb/in 2 occurred at the leading edge region 
(figs. 13 and 14). The orbiter wing skin buckling stresses are approximately a x * 

-12,000 lb/in 2 and Oy * -25,000 lb/in 2 . The peak value of the shear stress J T x | * 

7877 lb/in 2 was located at the wing root trailing edge zone (fig. 15). Figures 16 
to 18 respectively show the distributions of o x , Oy, and x X y in the orbiter wing 

upper skin. The peak tension of a x = 3087 lb/in 2 occurred near the wheel well, and 
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the peak compression of a x - -2947 lb/in 2 was located at the wing root trailing edge 

region (fig. 16). For Oy, the peak tension 2483 lb/in 2 was located at the midspan 

trailing edge region, and the peak compression ( Oy « -1371 in/lb 2 ) occurred at the 
first bay of midspan (fig. 17). 

Like the orbiter wing lower skin, the maximum shear | i X y | » 2137 lb/in 2 was 
located at the wing root trailing edge region. Figure 19 shows the shear stress 
distribution in the wing spars, leading edge panels, wing root wall, and the ele- 
von support panels. The peak shear occurred at the trailing edge region of the 
wing root wall (fuselage wall). 

Figure 20 shows the axial stresses in the orbiter wing spar caps, wing root rib 
caps, and other rod elements. The maximum axial tension of 7988 lb/ in 2 occurred in 

the aft wheel well wall lower spar cap; the peak axial compression of -15,408 lb/ in 2 
occurred at the wing root rib lower cap near the midchord region. Figure 21 gives 
the axial stresses in the orbiter wing rib truss members. The peak tension of 
5412 lb/in 2 occurred at the vertical truss member of the last bay of the rib next to 

to the wing root rib; the peak compression of -9038 lb/in 2 occurred at the lower 
horizontal truss member of the wing root rib. 

Figures 22 to 24 respectively show the chordwise distributions of three thermal 
stresses o x , Oy, and T X y at wing cross section Yq- 240. These stresses were calcu- 
lated using the WS240 and WING structural models. For the WS240 model, two sets of 
curves were presented. The first set (solid curves) is for the "plane stress" bound- 
ary condition and the second set (dotted curves) is for the "plane strain" deforma- 
tion. In figure 23, the distribution of Oy calculated from WS240 model for the wing 

lower skin and upper skin of bay 1 exhibits stress release zones near the center 
region of each bay. This implies that the wing skins in those regions .have bulged 
out because of thermal loading. The WING model did not have the capability to show 
the above stress release effect because of insufficient finite elements. Poor corre- 
lation between the predictions based on WS240 and predictions based on WING models 
could be attributed to: 

/ 

1. The two.sets of boundary conditions used for the WS240 model may 
not represent the actual deformation conditions. 

2 . The elements used for the WING model could be too coarse to give 
more accurate stress distributions. 

The variations between the calculated stress data of figures 22 to 24 illustrate 
how sensitive the wing segment models are to boundary conditions. More importantly, 
the differences between the calculated stresses from the wing segment model and the 
wing model are significant in terms of percentage and sign. Even though the wing 
model is relatively coarse, it is still considered to be the most accurate anal- 
ysis (of those presented here) because it considers temperature and stress effects 
throughout the entire wing. The WING model is also believed to be the least sus- 
ceptible to those boundary conditions imposed in the case of the wing segment model. 
An increased number of elements in the entire wing model should improve the distri- 
bution of calculated stresses between spar caps, but the magnitude and range of 
stresses would be expected to change very little. 
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Figure 25 shows the distribution of the axial stress a x in the fuselage struc- 
ture. The peak compression of a x = -7653 lb/in 2 in the fuselage bottom skin occurred 

near the vertical wall, and the peak tension of a x = 7385 lb/in 2 occurred near the 

glove leading edge region. Tables 1 and 2 summarize the major peak stresses that 
occurred in the orbiter major structural components. 

Figure 26 shows the deflection curves of the orbiter wing leading and trailing 
edges. The wingtip deflection induced by the thermal loading is 0.94 in. The inner 
span region of the wing is twisted 'in the direction of increasing angle of attack, but 
the outer span region of the wing is twisted in the opposite direction. Figure 27 

shows the spanwise change in angle of attack Aa as a result of the thermal loading. 

The maximum value of j Aa j * 0.1° occurred near the wingtip. 

Figure 28 shows the deformed shapes of the wing cross section at Yg-240 predicted 

by using the wing segment model under the two boundary conditions mentioned earlier. 
The entire wing cross section has curved up slightly because of thermal loading. The 
deformation at bay 1 is most conspicuous for boundary condition 1 (fig. 28(a)), with 
a vertical displacement of 0.110 in. at the bay 1 lower skin (honeycomb-core sandwich 
skin). For boundary condition 2 (fig. 28(b)), the peak vertical displacement of 
0.121 in. occurred at the bay 1 upper skin. Figure 29 shows the deformed fuselage 
cross section at Xq 877. Tne thermal loading tended to flatten the fuselage outer 

bottom, and caused the glove leading edge to move slightly upward, with upward dis- 

placement of 0.013 in. 


CONCLUS IONS 


* Thermal stress analyses were performed on the space shuttle orbiter subjected to 
STS -5 reentry thermal loading. The whole wing, one midspan wing segment, and one 
midfuselage segment were selected for the analyses. The whole wing was found to be 
twisted under the thermal loading. The inboard span region of the wing is twisted 
in the direction of increasing angle of attack, but the outer span region of the wing 
is twisted in the opposite direction. The maximum twist angle of the wing is about 
0.12°. The correlation between the stress predictions using the whole wing model and 
the wing segment model was rather poor. The reasons for this may be (1) the boundary 
conditions used in the wing segment model may not accurately represent the actual 
deformation field in the whole wing and (2) the whole wing model could be too coarse 
to give sufficiently accurate stress distributions. Even though the wing model is 
relatively coarse, it is still considered to be the most accurate analysis compared 
with the wing segment model because it considers temperature and stress effects 
throughout the entire wing. The fuselage cross section was found to deform in such 
a way that the fuselage outer bottom was slightly flattened and the glove leading 
edge moved slightly upward. Finally, it was found that the thermal stress levels 
induced in the orbiter structure were relatively low. 
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TABLE 1. - PEAK STRESSES IN ORBITER STRUCTURAL COMPONENTS 




lb/in 2 

(7y t 

lb/in 2 

| i xy | , lb/in 2 

Structural 

component 

Peak 

tension 

Peak 

compression 

Peak 

tension 

Peak 

compression 

Peak 

shear 

Wing lower skin 
(Fig. 13 to 15) 

2836 

-9097 

2405 

-2897 

7877 

Wing upper skin 
(fig. 16 to 18) 

3087 

-2947 

2483 

-1371 

2137 

SPAR web/fuselage 
wall (fig. 19) 

— 


— 


1829 

Fuselage skin 
(fig. 27) 

7411 

-7628 

533 

-574 

69 


TABLE 2. - PEAK AXIAL STRESSES IN ORBITER 
SLENDER STRUCTURAL COMPONENTS 


Structural 

component 

Axial tension, 
lb/in^ 

Axial compression, 
lb/in 2 

Wing spar caps and 
wing root rib caps 
(fig. 20) 

7988 

-15,408 

Wing rib truss 

members (fig. 21) 

5412 

-9,038 



Figure I. Regions of space shuttle orbiter modeled . 
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Figure 6 . Space shuttle orbiter wing SPAR structural model. TPS , wheel 
well door , and landing gear excluded • 




Y.-W 



WS S2-VO 


2o4 a«io points 
I 2. 1 CQU4D2 
139 Ei_iar^r*MTS 

Figure 7. NASTRAN structural iaodel for orbiter wing segment WS240 . 
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Figure 8. NASTRAN structural model for orbiter fuselage segment FS877 , 


6 







TEMPERATURE. °C TEMPERATURE 


O STS-5 FLIGHT DATA 

SPAR, 100% TPS THICKNESS 

TOUCHDOWN TIME, 1821 sec 



Figure 10. Orbiter fuselage skin temperature — time histories at station 
Xq 877 . STS-5 flight (ref. 5). 
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Figure 11. Chordwise distribution of aluminum skin temperatures at 
uring station | Yq j 240. Time * 1700 sec. STS-5 flight. 
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Figure 13. Distribution of normal stress o x (in lb/in^J in the orbiter 
wing lower skin. STS-5 thermal loading. Time * 1700 sec. 
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Figure 14. Distribution of normal stress Oy (in lb/in? ) in the orbiter 
wing lower skin. STS-5 thermal loading. Time m 1700 sec. 
































































Figure 15 • Distribution of shear stress i X y (in Ib/in 2) in the orbiter 
wing lower sJcin* 5T5-5 thermal loading • Time * 2700 sec* 
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Figure 17. Distribution of normal stress o y (in lb/in?) in the orbiter 
wing upper skin. STS-5 thermal loading.. Time m 1700 sec. 
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Figure 18. Distribution of shear stress r X y ( in lb/ in 2 ) in the orbiter 
wing upper skin. STS- 5 thermal loading . Time * 1700 sec. 
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Figure 20. Axial stresses (in lb/in 2 ) in the orbiter wing spar caps, wing 
root rib caps, and other rod elements. STS-5 thermal loading. Time m 
1700 sec . 
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